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Evaluation of the Exponential Integral for Large Complex 

Arguments 

John Todd 

Two methods of evaluating the exponential integral for (large) complex arguments are 
discussed. It is shown that the Laguerre quadrature method is more efficient than the 
asymptotic expansion. Examples are given to show the practicability of the Laguerre 
method. 

1. Introduction 
An extensive table of the exponential integral 



Ei(z)= I (e~ u /u)du z = x + iy 



has been prepared by the National Bureau of Standards [1]; l the introduction to the table 
gives a precise definition of this function. This table covers the range |x|<20, |?/|<20, with 
arguments variously spaced. In order to compute Ei(z) outside 4 this range, (or within this 
range, at points where interpolation is awkward, which are those with comparatively large 
argument), several methods can be considered. We examine two in detail and recommend the 
Laguerre quadrature, which is certainly practicable in the case of isolated arguments. 

2. Asymptotic Series 

The first of these is the use of the asymptotic series (see e. g., [1]), 

wt N _ z (\ 1! . 2! 3! . \ 

The remainder after n terms in the expansion for e z E x (z) can be obtained by integration by parts 
and is 

If we put u — z-\-p, where p is real, then 

y ' Jo [(x + P )+iy] n+1 
which can be estimated as 

\B„\<S n =n\[x 2 +y 2 ]- iin+1) if «>0 

\B n \<T n =n\\y\- n - 1 if x<0, 

from which the asymptotic character of the series is evident. In fact, 

S n+1 _ n + 1 
S„ (x'+y^' 

and so S n decreases as n increases until n exceeds (x 2 -{-y 2 )K To obtain the least the least value 
of S n we have to take n=(x 2j [-y 2 )*=\z\. For this value of n, we find, using Stirling's formula, 
that S n is about (2w/n)h~ n . In the second case, the least value of T n is formally the same, but 

1 Figures in brackets indicate the literature references at the end of this paper. 
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now n=|i/|. The decrease in the effectiveness of this representation as z approaches the 
negative real axis is apparent from the estimates. 

We note that the evaluation of the sums of the first n terms of the series for e z Ei(z), as a 
polynomial in z' 1 , by successive multiplication by z~ l for instance, involves about An real 
multiplications, z being complex. 

3. Laguerre Quadrature 

The second method is the use of a Laguerre quadrature. We shall show that this is more 
efficient than the asymptotic series method just described. 

We have 

eFExiz) = Ii-U 2 , 
where 



f ro x + p J T C V J 

Jo (x + p) 2 +y 2 p ' Jo (x + p) 2 +y 2 

It is well known [3] that, for any n, 



p. 



■/; 



where the xl n) are the zeros of the Laguerre polynomial L n (t), and the X/ n) are the correspond- 
ing Christoffel numbers. The zj* y , \; (n) have been tabulated for n=l(l)15, i=l(\)n by 
Salzer and Zucker [2]. 

Moreover, it is known [3] that for some £, < £ < °° y 

We shall now estimate this error. It is convenient to handle the two integrals separately. 
We calculate the 2nth derivatives of 

and find 

2? l = (7i!)«rr a »- 1 cos (2n+l)0 u 

R 2 ^(n\) 2 r2 2 *- 1 sin (2n+l)$ 2 , 

where the £*=£*(#, y) satisfy 0<{<< » f r i =[(x J r ^ i ) 2 +y 2 ]^, and Z - is defined by cos Q i ={x+% l )lr i) 

sin 6 i =y/r i . 

Hence 

1*1 ^i^ »*o. 

I^|<||^ *<o. 

The difficulties that were noticed in section 2, as z approaches the negative real axis, are still 
present. As before, there is an optimum value of n, which is about \z\ in the first case, \y\ 
in the second. With this choice of n, the value of our bounds for the R t are about 27re~ 2w in 
each case. 

The important fact to notice is that our estimate for the error committed by using an 
7i-point Laguerre quadrature is approximately the square of that for the estimate of the error 
committed by neglecting all terms in the asymptotic series after the nth. This, together with 
the fact that about the same number (An) of real multiplications are required in both methods, 
indicates that the Laguerre method is more efficient. 

314 



4. Examples 

In order to demonstrate the practicability of the Laguerre method, and to show that our 
conclusions are not dependent on the crudity of our appraisals of the remainders, we give three 
examples in detail. We choose n=5, which is optimum in one but not in the others. For 
completeness, we repeat (from [2]) the values of x ( (5) , \* (5) . 



x} 5) 


X, <5) 


0. 26356 03197 18 


0. 52175 56105 83 


1. 41340 30591 07 


.39866 68110 83 


3. 59642 57710 41 


. 07594 24496 817 


7. 08581 00058 59 


.00361 17586 7992 


12. 64080 08442 76 


.00002 33699 723858 



The superscripts (5) will be dropped from now on. 

We shall work to about 10 decimals throughout our examples. We begin with a case of 
z real: 2=10 and follow with £=l + 10i, and then z— — 10 + 5i; it is only in the last case that 
we are, with n=5, near the optimum number of terms. 

(a) 2=10. 

In this case/ 2 = 0, and we have only to deal with 



H>'-(iw^ 2x *-iofe =<21 - 



A complete worksheet for this calculation follows. The sum Qi is cumulated on the machine, 
there being no reason to record individual products. 

0^ = 10+2* a' 1 

10. 26356 03197 0. 09743 20771 

11.41340 30591 .08761 62872 

13. 59642 57710 . 07354 87412 

17. 08581 00059 . 05852 81002 

22. 64080 08443 . 04416 80490 

Q 1 = J2\ i ar 1 = 0. 09156 33319 

The correct value of e^E^lO) is 0.09156 3334. The actual error is about 2X10" 9 , 
whereas our estimate is 144X10" 9 . Using five terms of the asymptotic series, we obtain the 
value 0.09164. There is thus an actual error of about 8X10" 5 , compared with an estimated 
error of 12X10" 5 . 

(b) z=l + l0i. 

In this case we have 

e 1 + m E 1 (l + 10i)=I 1 -U 2 , 
where 

7l= Jo e "''(i+o 2 +ioo'* # ^ x< \i+^ 2 +^oo =S(1+Xi) ((i+^) 2 +ioo| 



and 



Jo '"'(1 +0 2 +100*^^ SX< '(l-fx 4 ) 2 +100 _Sl0 l(l+x < ) 2< +100r 
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A complete work sheet for this calculation follows. It is convenient to compute the common 
factors ji, in braces { } above, and then obtain the sums S a *y* an( i SIOt* by accumulation on 
the machine. 

fi^lOO + a! 

101.59658 47 

105.82451 43 

121. 12712 99 

165.38032 35 

286.07144 76 



a t =l+Xi 
1.26356 0320 
2. 41340 3059 
4.59642 5771 
8.08581 0006 
13. 64080 0844 



0. 00513 55625 
.00376 72444 
.00062 69648 
.00002 18391 
.00000 00817 



^ 1+lw -Ei(l + 10i)==Sam=0.01864 0471 
^ 1+lw JBi(l + 10i)=-Sl07<= -0.09551 6925 

The correct value of e 1+m Ei(l + 10i) is 0.01864 049-0.09551 688^. 

The actual error is of the order of 5X10 -8 , whereas our estimate is about 144X10" 9 . 

five terms of the asymptotic series gives 

0.018595-0.095600i 



Using 



There is thus an actual error of about 1X10 4 , compared with an estimated error of about 
1.2X10" 4 . 



(c) 2= — 10 + 5i 



a i = — 10 + x i 
-9.73643 9680 
-8. 58659 6941 
-6.40357 4229 
-2.91418 9994 
2. 64080 0844 



/3<=25+a? 
119.79825 76 
98.72964 70 
66. 00576 29 
33.49250 33 
31.97382 91 



0.00435 52855 
.00403 79645 
.00115 05428 
.00010 78378 
. 00000 07309 



#«- 10+M #i(-10 + 5i)=2><7< = -0.08475 7264 

y«- 10+M Si(-10 + 5i)= — S57*= -0.04826 1807 

The correct value of e- 10+M £i(— 10 + 5i) is 

-.08475 749-0.04826 039i. 

The actual error is of the order of 2X10" 6 , whereas our estimate is about 3X10" 4 . 
five terms of the asymptotic series, we obtain the value 

-0.084855-0.048277;. 



Using 



There is thus an actual error of about 1X10 -4 , compared with an estimated error of about 
8X10" 3 . 

5. Remarks 

1 . A third method of attack is conceivable : The integration in the complex plane of the 
differential equation for Ei(z), taking z~ l as the independent variable, extending the ideas of 
Fox and Miller [4]. 

2. It is conceivable that the Laguerre method would be effective in cases when the integrand 
is not precisely of the form e~ l j(t). A case of importance might be that of the generalized 

exponential integral I (e~ w /w)du, where w=(a 2 +u 2 ) i . In such a case it would be conven- 
ient to follow the device of A. Reiz [5] and use the quadrature in the form 



J"no*=SM« w iw). 
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where /u/ n) =X/ tt) exp x$ n) ; the jui n) have been tabulated by Salzer and Zucker [2]. 

3. The table of Mashiko [6] gives E x (z), or rather auxiliary functions from which this can 
be obtained, for z=re ie when f-^O.OltO.OljO^ and 0=O(2 O )6O°(1 O )9O°. 
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